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Abstract 

Lemaitre-Tolman-Bondi (LTB) solutions are used frequently to describe the 
collapse or expansion of spherically symmetric inhomogeneous mass distributions 
in the Universe. These exact solutions are obtained in the synchronous gauge 
where nonlinear dynamics (with respect to the FLRW background) induce large 
deviations from the FLRW metric. In this paper we show explicitly that this 
is a gauge artefact (for realistic sub- horizon inhomogeneities). We write down 
the nonlinear gauge transformation from synchronous to Newtonian gauge for a 
general LTB solution using the fact that the peculiar velocities are small. In the 
latter gauge we recover the solution in the form of a weakly perturbed FLRW 
metric that is assumed in standard cosmology. Furthermore we show how to 
obtain the LTB solutions directly in Newtonian gauge and illustrate how the 
Newtonian approximation remains valid in the nonlinear regime where cosmo- 
logical perturbation theory breaks down. Finally we discuss the implications of 
our results for the backreaction scenario. 
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1 Introduction 



The Friedmann-Lemaitre-Roberston- Walker (FLRW) framework is the cornerstone of 
modern cosmology. Its key assumption is that our large scale Universe is isotropic and 
homogeneous up to small perturbations. Until not so long ago this was really an as- 
sumption, but fortunately, present observations like the Sloan Digital Sky Survey show 
that for distance scales > lOOMpc one does indeed seem to find statistical homogeneity 
At the same time however we know that at small scales < lOMpc our Universe 
looks like anything but the idealized FLRW Universe. Nearly all of the matter has 
clumped into clusters and galaxies with local density contrasts and dynamics that are 
clearly beyond a linear description on a FLRW background. Given the nonlinear na- 
ture of GR, one may wonder then if the FLRW description of our large scale Universe 
is really justified. And if the corresponding scale factor evolves according to the Fried- 
mann equations, with an effective matter density that is simply the averaged actual 
matter density. 

This issue was studied in the past but gained much interest over the last few years 
with the discovery of the apparent cosmic acceleration (see [21 E] for a review and an 
extensive list of references). Within the conventional FLRW framework this accelera- 
tion demands the presence of a cosmological constant, some other form of dark energy 
or a large distance modification of gravity. All cases require the ad hoc introduction 
of a new mass scale in the Lagrangian that is suspiciously of the same magnitude as 
the present Hubble constant. In light of this, several papers advocate the so called 
backreaction scenario [3J, HJ El El El El E] ■ According to this scenario the conventional 
FLRW framework is not justified. The hope is then that a correct averaging of the 
inhomogeneities at small distance scales would lead to an effective energy-momentum 
tensor to be used in the FLRW description at large scales, that would obviate the need 
for dark energy or modified gravity. A nice feature of this scenario is that it would 
naturally solve the coincidence problem (why now?) since it ties the epoch of cosmic 
acceleration to the epoch of nonlinear structure formation. 

However, a strong argument against this scenario was raised by Ishibashi and Wald. 
In [10] they simply point out that, despite the nonlinear dynamics at small scales, the 
actual space-time metric of our Universe seems to be very well approximated by a 

1 That is if we don't consider the possibility of us living at the center of a Hubble-sized void. 
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Newtonianly perturbed FLRW metric at all scales. By this, they mean a metric of the 
formal 

ds 2 = - (1 + 2ip) dt 2 + a(t) 2 (1 - 2^0 (dx 2 + dy 2 + dz 2 ) , (1.1) 

with 

H « i , {d t ipf « \d4de\> , (d^) 2 « (didjWdidrf) (1.2) 

or 

If (11.11) is really a good approximation everywhere to the metric of the Universe, 
then it is evident that the backreaction scenario does not work. Indeed, in that case 
one can immediately show that the scale factor a(t) in the metric evolves according to 
the usual Friedmann equations, up to small corrections. Of course what we observe 
are rays of light, not a metric. But given a metric of the form (II. ip one can resort to 
a standard lensing analysis to show that the effect of the small scale inhomogeneities 
on cosmological observations is indeed small and can not explain away dark energy. 
(See [TT] for instance, for an estimation of the effect in the context of supernovae 
observations.) 

Now, one thing is the observation that the nonlinear dynamics at small scales do 
not preclude a weakly perturbed FLRW metric. Another, stronger statement, is that 
the metric (11.11) actually is a good approximation to the exact metric of the Universe. 
To argue the latter one has to resort to an expansion in the peculiar velocity v of the 
full set of Einstein equations. In this expansion, that lies at the root of the Newtonian 
approximation in comoving coordinates, one finds that (11.11) is indeed the leading order 
approximation to the metric and that the higher order terms are small. However, 
this perturbative argument could in principle still be invalidated by non-perturbative 
effects. 

In this paper we study this issue in the case of spherically symmetric space-times, 
the advantage being that we can use the exact LTB solutions [12J to check the validity 
of the Newtonian approximation. Not so surprising, the conclusion will be that the 
Newtonian approximation is valid both in the linear and nonlinear regime. As such 
this serves as a nice illustration of the compatibility of a global Newtonianly perturbed 
FLRW metric with local nonlinear dynamics. A subtlety will be that the LTB solutions 
are formulated in the synchronous gauge. In this gauge, starting from small initial 
2 For simplicity we are assuming flat FLRW here. 
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perturbations, the weak field description of the metric breaks down at the same time 
when the dynamics go nonlinear. In the next section we will show this explicitly for the 
LTB solutions, by linearizing the exact solutions on the FLRW background to recover 
the standard results of cosmological perturbation theory. In section 3, we will show 
however that this phenomenon is specific to the synchronous gauge. We will explicitly 
demonstrate that, for small peculiar velocities, the corresponding metric in Newtonian 
gauge remains of the weakly perturbed FLRW form throughout the entire evolution. 
In doing so, we will also verify that the solution that one obtains in the Newtonian 
gauge by solving the equations in the Newtonian approximation, agrees very well with 
the corresponding exact LTB solution. 

The gauge transformations from synchronous to Newtonian gauge have been con- 
sidered before in the context of LTB solutions. A valid application of the linear trans- 
formations was used in [13], an invalid application can be found in [14] . In [15] the 
full nonlinear transformation was performed on a specific LTB solution, obtaining a 
different result than ours. We will comment on this further on. 

2 LTB solutions 

The LTB solutions are exact solutions to the Einstein equations for spherically sym- 
metric distributions of irrotational pressureless dust, so with energy- momentum tensor 
T[j,v(r,t) = pUpU v . They are obtained in the synchronous gauge (goo = — l,<?oi = 0) 
and from now on we will label the radial and time coordinate in this gauge as (f, t), 
reserving (r, t) for the Newtonian gauge. Notice also that for pressureless dust, the 
synchronous gauge coincides with the comoving gauge (m m = (1, 0, 0, 0)). Explicitly the 
solutions look like B 

ds * = _ d p + Rl2 ( r 'f dr 2 + R(r } tfdn 2 , (2.3) 
1 + 2E{r) 

3 We use units such that c = 1 . 
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witlfi 



«!M +fl(r - il% - (1 = «?M_w) = m (24) 

2 v ; v ; 2 R(r,t) K ! V ' 

, ~, M'(f) , \ 

= *mim- (2 - 5) 

solved by (with E(r) G R): 

i2 = (cosh (wv 7 ^) - l) , (2.6) 

GjvM .sinh(wV2~B) N ^ „ N 

A particular solution is then specified by some set of initial conditions 
{p(f, tj), i?(f , tj), R(r, ti)} at time that translate to a certain choice for the free func- 
tions {.E(f), M(f), £{,(f)} via eqs. (12 .41) - (12. 7p . One of the initial conditions fixes the 
residual gauge freedom, we will choose R(r,ti) = Ojf , where a« = a(tj) is the scale 
factor at time ti of the FLRW background at infinity. The two other initial conditions 
fix the initial density and peculiar velocity v, where the latter is defined as (H = a /a): 

v = R-HR. (2.8) 

Some comments are in order before moving on. For the remainder of the paper we 
will only consider solutions that approach the Einstein- de Sitter solution (flat, matter 
dominated FLRW, with a(t) = (t/to) 2 ^ 3 ) at infinity, but it's straightforward to extrap- 
olate our results to situations with non-flat FLRW backgrounds. One could even add 
radiation and/or dark energy for the FLRW background at infinity by taking slightly 
modified LTB solutions. Furthermore we will restrict ourselves to solutions that at 
some early time were close to the FLRW solution, in accord with the standard picture 
of the early Universe. It's also clear that v(r, t) corresponds to the velocity of the shell 
of matter labelled by f with respect to the FLRW background. Of course, by definition, 
the peculiar velocity in comoving/ synchronous gauge is zero. But we will show later 
on that v does indeed correspond to the proper peculiar velocity in Newtonian gauge 
(at least if v <C 1). 



4 Some conventions for the derivatives: the prime stands for a radial and the dot for a time derivative 
with respect to the arguments of the function at hand. Later on we will also use the coordinates (r, t) 
in Newtonian gauge, so we will have for instance a' = a'(r, t) = ^ or d r R = d r R(r, t) = Rd r t+R'd r r . 
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Let us now consider some particular solution, describing the evolution of an initial 
inhomogeneity, characterized by a velocity and density profile at some early time ti. 

Vi (f) = v(f,ti), (2.9) 
Pi(r) = p(r,ti) = Pflrw(U)(1 + 5(f,ii)) = p F LRw{ti){l + Si(f)) , (2.10) 

with 



5i,t)i<l , limt>i(r)=0 , \imvi(r)=0 , / drr z 5i(r) = 0. (2.11) 

f->0 f-»oo J Q 

From eqs. (I2.4p - fl2.7l) and the first Friedmann equation we immediately find: 

M{f) = —^p FLRW (ti)(l + 3I 5i ), (2.12) 

E{f) = Vi H t a l f-^(H l a i f) 2 I 5i + ^-, (2.13) 



E 

Hi V 5 idiHif) 



h{f) « —[- J -—- + I Si ) , (2.14) 



where we have defined 

I 5i (r) = ^ I dr'r%(r'). (2.15) 



It is instructive now to expand the exact solution for R(f,t) in 5i,Vi (treating both 
quantities as small parameters of the same order). This is precisely the expansion that 
is considered in cosmological perturbation theory to obtain an approximate solution to 
the equations of motion for more generic (non-symmetric) initial conditions. For the 
case at hand we find at leading order (keeping the full expression of E ~ Vi ~ 5i): 

t\ 2/S [2 E \ I, [2 E 



R( f , t) * M t) | 1+ | [jj ^ ^ ) +h ,- l{ ^— + h ^ + ... 

(2.16) 

with the dominant part (for t^>tj) of the higher order terms going like 

t\ 2/3 E , 

(2.17) 



ti) (cnHir) 2 

We recognize the first term in the linear term corresponding to the growing mode, 
and the last term corresponding to the decaying mode of cosmological perturbation 
theory (see for instance [16] ). Furthermore, from the behavior of the higher order 
terms we see that the expansion breaks down at a time t n i when 

U \ 2/3 c 



Now, as we already commented in the introduction, one thing is the breakdown of 
the expansion in Si, another thing is the breakdown of the metric description in terms 
of a weakly perturbed FLRW metric. But as one can see from f 1 2 . X 6 [) . the breakdown 
of the former automatically implies the breakdown of the latter for the synchronous 
gauge solution that we are considering. Indeed we find 

(2.19) 

a(t)r 

at t ~ t n \. In addition one can see from the exact expression for the matter density 
(12 .5|) . that the density contrast also becomes large at the same time: S(r,t n i) ~ 1. 

Notice that the converse is not necessarily true, one can have solutions with shell- 
crossing singularities (g?f — R' — 0), for which the description of the metric in terms 
of a weakly perturbed FLRW clearly breaks down, that are perfectly well described by 
the perturbative expansion in S{ for times arbitrary close to the time of shell-crossing. 
The "small u expansion" of [T7j uses precisely this feature in combination with the 
exact expression (12.51) to describe the density. 

3 From synchronous to Newtonian gauge (and back) 
3.1 The linear case 

We saw in the previous section how the description of the LTB metric as a weakly 
perturbed FLRW metric in synchronous gauge, breaks down at the same time when 
the dynamics on the FLRW background become nonlinear. Let us now see what 
happens in Newtonian (or Poisson) gauge. In spherically symmetric situations, given 
the standard angular coordinates, this gauge is specified by the conditions g rt = 
and g rr r 2 = ggg. So we want to find the coordinate transformation (f, t) — > (r,t) that 
changes the metric to: 

ds 2 = -dt 2 (1 + 2(j)(r, t)) + a 2 {t) (dr 2 + r 2 dQ 2 ) (1 - 2^(r, t)) . (3.20) 

At times t <C t n j, when the dynamics are still linear, the coordinate transformation 
will also be lineajfl and one can resort directly to the machinery developed in the context 
of cosmological perturbation theory to obtain the potentials [16J. This was done for 

5 By a linear coordinate transformation we mean that t = t + a and f = r + (3 with a', a, (3' ,(3 <C 1. 
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instance in [13] and for our conventions we find in a completely similar fashion that 
the resulting potentials read: 

0«V + W) 2 «-§ Jdr^l + ^(y) 5/ Y dr Q^ + r(a ^ )2 ^ (r) ) • (3 ' 21) 

One can again recognize the growing mode and the decaying mode. As we should, 
we find the former to result in a constant term for the potential in Newtonian gauge. 
This, in contrast with the situation for synchronous gauge, implies that the metric 
keeps its form of a weakly perturbed FLRW ~ 1) up to the time t n i, when the 
dynamics become nonlinear. Of course one can not resort to cosmological perturbation 
theory anymore to argue for the smallness of the fields at later times, when the fields 
and the actual density contrast 5, become nonlinear functions of S^. Still, from the 
argument of Ishibashi and Wald, we might expect the potential (ip ~ <p) to remain 
small for sub-horizon inhomogeneities. As we will show explicitly in the next subsection, 
this is indeed the case for the LTB solutions in Newtonian gauge, provided that the 
peculiar velocities remain small. For the LTB solutions this will hold as long as we stop 
the evolution well before a central singularity develops. In reality, when we consider 
the collapse of a cluster for instance, the (effective) pressure that arises during the 
virialization, will halt the collapse. 

3.2 The nonlinear case 

We could proceed now in the same way as we did for the linear case. So, starting 
from the LTB solution we could look for the coordinate transformation (f,t) — > (r,t), 
which is now nonlinear, to obtain the metric in the Newtonian gauge. This strategy 
was used in |15j . But we find it more instructive go the other way. We will start from 
the equations in Newtonian gauge for spherically symmetric situations, with metric 
( 13.201) . showing explicitly how the expansion in the peculiar velocity v gives rise to 
the familiar equations of the Newtonian approximation to GR. Then we will perform 
the coordinate transformation to synchronous gauge, demonstrating explicitly that the 
solutions obtained in the Newtonian approximation are indeed good approximations 
to the exact LTB solutions. 

First we should explain at last what we mean exactly by the expansion in v. From 
Newtonian physics one gets the following order of magnitude estimates in the case of 
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a non-virialized system with density contrast S and distance scale L[18j: 

2 

ij) ~ <p ~ 5(FL) 2 , 0/ ~ 5# 2 L , w ~ &HX , , t> ~ ^- , . . . . (3.22) 

When assigning a power of v to a particular term in the expansion, we will use the order 
of magnitude estimates above, with 8 ~ 1. This ensures that the expansion remains 
valid in the nonlinear regime 5 > 1. 

Let us now apply this expansion to the full set of Einstein equations. In the spher- 
ically symmetric case that we are considering, there are four independent equations: 

G tt ~ StiGmP i G tr ~ STiG^pva , G rr ~ 8nGNpv 2 a 2 , G w = sin 2 9Gee = , 

(3.23) 

for the perfect dust energy-momentum tensor = pu^u^, with ~ (l,v/a, 0,0). 
We easily get the fields and i/j from the (#, 9) and (t, t) equations. The former 

Gee = r V - ^ + V - V) + 0(v 4 ) = , (3.24) 

immediately tells us that ip = <p + 0(t> 4 )E The latter then reduces to the Poisson 
equation (using the FLRW equation for the background), with the potential sourced 
by the density contrast: 

V 2 <p = 47cG N a 2 p FLRW 5 +H 2 0{v 2 ). (3.25) 

The metric (13.201) will then solve the other two Einstein equations if the energy- 
momentum tensor is conserved. That is if p and v obey the familiar ideal fluid equations 
in comoving coordinates (for zero pressure) |18j : 

V M T^ =(p + 3Hp + -(pv)' + -pv) + -^-<D(v 3 ) = , (3.26) 
\ a ra J G^L 

* ( V T »r _ v _ v T »t\ = L + — + HV + 4-) +\o{v A ) = 0. (3.27) 
p V a / \ a a J L 

The solution in the Newtonian approximation can now be obtained by solving 
the eqs. (13.2 51) - 03.271) at leading order, for a particular initial velocity and density 
profile Vi(r),5i(r) <C 1 at some early time t{. One can verify that for these solutions, 
the order of magnitude estimates (I3.22p are correct, so the higher order terms in the 

6 We are imposing the boundary conditions <fi, -0, 4>' , -0' — ■> for r — > oo, to match the FLRW 
background at infinity. 
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expansion will indeed be suppressed, both in the linear (5 <C 1) and nonlinear (5 > 1) 
regime as long as v remains small. This type of argument on the validity of the 
Newtonian approximation could be criticized for being circular. In a sense we are 
using the Newtonian approximation to justify itself. But as we will now show by going 
to the synchronous gauge, the approximate solution obtained in the Newtonian gauge, 
is indeed a very good approximation to the exact LTB solution. 
We define the coordinate transformation by: 

R(f,t) 2 = a(f)V(l -2^(r,t)), (3.28) 
t = t + a{r,t), (3.29) 

for some functions R(r,t) and a(r, t) that will be determined from the synchronous 
gauge conditions g~ tf = and gn = —1. We will do this in the same expansion in v 
that we used to obtain the Newtonian approximation, retaining only the terms up to 
0(v 2 ). Immediately we can anticipate that 

a'tt-av, (3.30) 

from the condition that u r should be zero in the synchronous/comoving gauge. Apply- 
ing the coordinate transformation on the metric, with the use of 

dt = dt (1 + a) + dra' , (3.31) 
R'dr « dt (or - i?) + dr (a(l - if)) - ar^p' - Ra'^j , (3.32) 

we then find 

a'R' 

g if w + R (R - dr) = , (3.33) 

d 

for 

a' 

R-ar^R-HR^ (3.34) 

a 

Employing this expression it is then straightforward to show that g% w — 1, if v obeys 
the Euler equation (13.271) . So the transformation fl3.28p - fl3.29p with a obeying (13.301) 
and R obeying (13.341) indeed takes us to the synchronous gauge. 

We will now show that R approximately solves the LTB equation if v, p and <f) 
solve the eqs. fl3.25p - fl3.27p . that we obtained in the Newtonian approximation. Let us 
first look at the continuity equation (13.261) and show that it is equivalent to the LTB 
expression (12. 5p for the matter density. Expressing the (r, t) derivatives in terms of 
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(f, t) derivatives through eqs. (13.311) and (13.321) . this equation becomes (again keeping 
the appropriate powers of v): 

d iP + 3Hp + p% + — pv w . (3.35) 
R' ra 

Using the expressions (13.281) and (13.341) for R and v, this reduces to 

R 




aO, (3.36) 

solved by the LTB expression (12.51) for p. 

Employing this LTB expression for p in the Poisson equation (13.251) for <p, we find 
in a similar way that: 



r 2 V./ V R'R? 2 

aif 2 .R 



a 
R 2 



drG N M'(f 



GmM H 2 R 



If we now use this expression for 0' in combination with a(t) = (t/to) 2 ^ 3 , we finally find 
that the Euler equation (I3.27P reduces to: 

which is now indeed solved by the LTB equation (12.41) . So we have demonstrated 
that for a solution obtained in the Newtonian approximation one finds gee = R(f,t) 2 
in the synchronous gauge, with R an approximate solution of the exact LTB equa- 
tions. Keeping track of the omitted terms in our expansion one can show that the 
approximation holds up to terms ~ iA As for the initial conditions, it is easy to see 
that the initial density and velocity profile at time ti in the Newtonian gauge, trans- 
late to (approximately) the same initial conditions at time ti for the LTB solutions: 
5i(r),Vi(r) pa 5i(r),Vi(f). That is if r ~ f around the time ti, which is true if we fix 
the residual gauge degree of freedom by the condition i?(f, tj) = ajf, as we did in the 
previous section. 

A nice cross check of our derivation is provided by the calculation of g?r- From the 
exact LTB solution we know that we should find 

R' 2 

9ff = T+2E ' (3 ' 39) 
10 



whereas from the coordinate transformation fl3.3ip - fl3.32p we find, 



a' 2 



g ff « R' 2 1 + 2 -a'r + ip'r - — . (3.40) 



a 2 



Using the expressions for i? (I3.28|) and a (13.341) this indeed reduces to 



i? 



1? 



■/2 



l + 2£' 

where on the last line we have used the LTB equation (12.41) and fi< 1. 



(3.41) 



To recapitulate, we have demonstrated that the coordinate transformation, implic- 
itly defined by: 



R(r,t) 2 = a(t) 2 r 2 (l-2^(r,t)), (3.42) 

POO 

i = t + a(t) / dr'v(r',t), (3.43) 

J r 

with 

R(r, t) - H(i)R(f, t) = v(r, t) , (3.44) 

takes the metric (13.201) in Newtonian gauge to synchronous gauge, provided v is small. 
In the synchronous gauge we recover the metric of the LTB form (12. 3p . with R and 
p approximately obeying the LTB equations (12. 4p and (12.51) in synchronous coordi- 
nates, if p, v,(j) ~ if} are solutions of the Poisson, continuity and Euler equations in the 
Newtonian coordinates. 

It's clear that we can use our results also in the other direction. Namely if we start 
from the exact LTB solution in synchronous gauge, the transformation (13 .421) - (13 .431) 
will take the metric to a Newtonianly perturbed FLRW metric with p, if), v solving the 
Poisson, continuity and Euler equation in Newtonian coordinates at leading order in 
v. We should stress (at the risk of being repetitive) that our analysis is valid both in 
the linear and nonlinear regime, as long as the peculiar velocities remain small. Notice 
also that the conditions on if) (jl.2p imposed by Ishibashi and Wald are automatically 
satisfied, again if v is small. 

As we mentioned, in [15] a similar analysis was performed for the specific case of 
zero initial velocity and an initial density profile describing a constant over-density 
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in the core surrounded by a finite region with a constant under-density. Surprisingly, 
the authors found a considerable difference for the fields ip and when the dynamics 
go nonlinear, seemingly in conflict with the Newtonian approximation and with our 
results. However, it turns out that the reason for the discrepancy lies in a calculational 
error [19]. 

4 Conclusions 

In this paper we have demonstrated explicitly how one can recover the exact LTB 
solutions in synchronous gauge, from the corresponding solutions obtained in the New- 
tonian approximation. This was done by applying the full nonlinear coordinate trans- 
formation from Newtonian to synchronous gauge, using the fact that the peculiar ve- 
locities remain small, which is the case for realistic sub-horizon inhomogeneities. 

As we have illustrated, in their original form the LTB solutions display a breakdown 
of the weak field description of the metric, when the dynamics go nonlinear with respect 
to the FLRW background. This has been used to argue for the backreaction scenario by 
Rasanen [3l 0] for instance. Indeed, once the weak field description breaks down, one 
would not a expect a priori that a Universe full of collapsing structures and expanding 
voids would still give a FLRW metric on the average with a normal evolution of the 
scale factor. However, as we have shown, the breakdown of the weak field description 
of the metric is in fact specific to the synchronous gauge. In Newtonian gauge, the 
LTB solutions are perfectly well described by a Newtonianly perturbed FLRW metric 
and it becomes straightforward to show that both the metric and the dynamics will on 
the average behave as a conventional FLRW Universe (up to small corrections). 

LTB solutions have been used also in the literature to model the effect of inho- 
mogeneities on light propagation in the so called Swiss cheese models [20J. Again, 
given the breakdown of the weak field description, one would not expect a priori to 
approximately recover the standard FLRW luminosity-distance redshift relation, when 
averaging over all directions. Yet this is precisely what is found and is again easy to 
understand from the weak field expansion in Newtonian gauge, as was demonstrated 
for LTB solutions in the linear regime in [13] . In a forthcoming publication we will use 
the full transformation (13 .421) - (13 .431) to analyze the Swiss cheese model of [21] in the 
nonlinear regime. We will find that the metric indeed reduces to a weakly perturbed 
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FLRW metric in the Newtonian gauge, which was used implicitly in the recent paper 
[22]. 

All this illustrates the strength of the Newtonian approximation in comoving co- 
ordinated in justifying the conventional FLRW framework. And it is clear that any 
serious backreaction scenario should explain where and how the Newtonian approxima- 
tion breaks down. Actually we know that the Newtonian approximation breaks down 
in the vicinity of black holes, and to our knowledge it has not been fully demonstrated 
yet that the corresponding backreaction is negligible. Notice also that in this paper 
we have only strictly proven the validity of the Newtonian approximation for spheri- 
cally symmetric space-times, that approach the FLRW solution at infinity. But as we 
commented in the introduction, for more general space-times one can show that the 
Newtonian approximation is at least self consistent. 

Finally we should comment on the local void scenario, which is another approach 
that involves the use of inhomogeneity in trying to dispose of dark energy. In this 
scenario one puts us near the center of a large void, typically described by an LTB 
solution. The mismatch between the local and global expansion can then explain the 
supernovae data [231 12H [25] . However, to explain other data sets, like those on the CMB 
and the large scale structure, one has to introduce additional features in the primordial 
power spectrum and the matter composition of the Universe [26l[27j . Another contrast 
with the backreaction scenario is that this scenario can be perfectly well described in 
a weak field description, with the local Hubble flow encoded in the Doppler term (see 
for instance [13]). Indeed, a good fit to the supernovae data requires a shift in the local 
Hubble parameter of order 10% with respect to the global value that is recovered at a 
distance r ~ lGpc. This translates to a maximal peculiar velocity v ~ AHr < 0.05 
which is still rather small and therefore validates the Newtonian approximation. This 
issue was studied by the authors of [H], for the type of void proposed in [24J. For 
this void, the present peculiar velocity is indeed small, and since the solution is in the 
nonlinear regime, one needs the full nonlinear transformation fl3.42p - fl3.43p to recover 
the metric as a Newtonian perturbation of the FLRW metric. Also, as was noticed in 
[T4] , a specific feature of this void solution is that the peculiar velocities become large 
(v > 1) in the past, at redshifts z > 500. At that time the void represents a nonlinear 



7 As opposed to the Newtonian approximation in physical coordinates, which breaks down at cos- 
mological distance scales. 
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super-horizon density fluctuation and the weak field description in Newtonian gauge 
will of course break down. 
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